Multiferroicity and spiral magnetism in FeV04 with quenched Fe orbital moments 
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FeV04 has been studied by heat capacity, magnetic susceptibihty, electric polarization and single 
crystal neutron diffraction experiments. The triclinic crystal structure is made of ^-shaped clusters 
of six Fe^+ ions, linked by V04~ groups. Two long-range magnetic ordering transitions occur at 
Tjvi=22K and Tjv2 = 15K. Both magnetic structures are incommensurate. That stable below Tjvi is 
coUinear with amplitude modulated moments whereas below Tjv2 the arrangement is non-coUinear 
with a helicoidal modulation. Below Tjv2, FeV04 becomes weakly ferroelectric coincidentally with 
the loss of the coUinearity of the magnetic structure. We conclude that FeV04 provides another 
example of frustrated spiral magnet similar to the classical TbMnOa compound. However, FeV04 
has quenched orbital moments and a particular structure clarifying the respective role of anisotropy 
and magnetic frustration in this type of multiferroic materials. 

PACS numbers: 25.40.Dn, 75.25.+Z, 77.80.-e 



There has been a recent surge of interest in a novel 
class of multiferroic materials in which ferroelectric- 
ity arises below a magnetic phase transition and as a 
direct consequence of complex magnetic ordering in sys- 
tems with strong magneto-electric inter actions [2[- This 
phenomenon, which requires lowering of the magneto- 
crystalline symmetry to a polar group, can only ap- 
pear in complex magnetic structures, typically stabilized 
by magnetic frustration or strongexchange competition, 
as in Ni3V2 08[i| and TbMnOsji]. By contrast, these 
"novel" or "type 11" multiferroics are simple binary 0| 
or ternary oxides — a desirable feature for prospective 
applications. 

In addition to fulfilling strict magneto-structural sym- 
metry requirements, type-II multiferroics must possess 
a microscopic mechanism to generate electric dipole mo- 
ments. A variety of such mechanisms have been proposed 
: magneto-striction is the only allowed mechanism for 
acentric collinear structures, and is active in the Ising sys- 
tem Ca3(Co,Mn)206 and, most likely, in the commen- 
surate phase of YMn2 05[i|. Most other "typc-II" multi- 
ferroics are cycloidal magnets, where non-coUinear spins 
arc key ingredients in the context of the so-called spin- 
current model [^, For this, a crucial role is played by 
rclativistic spin-orbit interaction, which can take place 
at the ligand ionic site, as for pure Cg systems, within 
the transition-metal <2g orbitals or between and eg 
orbitals Q . A particularly interesting case of the latter 
is provided by high-spin systems (S=5/2) where L=0 
in the free ion and the orbital angular momentum is sup- 
posedly absent. As for the acentric magnetic structure 
itself, this typically comes about when a connected net- 
work of superexchange (SE) interactions is destabilized 
by the presence of cither strong next-nearest neighbor in- 



teractions or geometrical frustration [J, [lH . 

In this letter, we describe a new multiferroic com- 
pound — FCVO4 — in which the magnetic ion is orbitally 
quenched Fe^+ (d^, L=0, S=5/2). The magnetic and 
dielectric phase diagram of FeV04, as determined from 
magnetization, specific heat and neutron diffraction mea- 
surements, is that of a typical cycloidal magnet: ferro- 
electricity appears below TAr2=15K, coinciding with the 
appearance of a non-coUinear incommensurate magnetic 
structure (phase II), whereas a second collinear incom- 
mensurate magnetic phase (I), stable between TAri= 22K 
and TAr2, is not ferroelectric. Uniquely, FeV04 docs not 
contain connected magnetic direct exchange or SE paths, 
and the magnetic modulation is primarily determined by 
a network of super-super-exchange (SSE) interactions. 
These pathways form loops connecting an odd number 
of Fe'^+, suggesting that frustration plays a key role in 
promoting non-coUinearity and ferroelectricity. 

Polycrystalline samples and single crystals were pre- 
pared following the procedures described in Refs IH and 
I13I . Electric measurements and specific heat and mag- 
netic susceptibility measurements were carried out on 
dense pellets of polycrystalline FeV04 using a PPMS 
Quantum Design cryostat. The electrical polarization 
was derived from integration of the pyroelectric current 
measured with a Keithley 6517 electrometer. Magneti- 
zation measurements on single crystals were made with 
a vertical-field SQUID magnetometer using a horizon- 
tal axis sample-rotator. Single crystal neutron nuclear 
and magnetic Bragg peak intensities were collected and 
at T=2K and T=18K on the four-circle diffractometer 
D15 (A = 1.1 74 A) at the Institut Laue-Langevin (ILL), 
France. All the nuclear and magnetic structure refine- 
ments were carried out with the program FullProFfl^. 
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FIG. 1: (Color) (a) Crystal structure and magnetic exchange 
frustration in FeV04. The basic magnetic unit are S-shaped 
clusters made of 6 Fe^"*" atoms (red, blue, and green large 
balls) weU spaced by tetrahedra (blue). The Fe-0 

bonds are drawn as thin grey lines (O atoms are not shown 
for clarity). The thick red lines represent intra-cluster inter- 
actions, with contributions from direct Fe-Fe exchange and 
Fe-O-Fe super-exchange. Iron atoms belonging to different 
clusters are coupled by Fe-O-O-Fe SSE paths mediated by 
the edges of the V04~ tetrahedra giving rise to effective inter- 
cluster Fe-Fe interactions (thin colored lines) 



The crystal structure of FeV04 is shown in Fig. [T] 
To facilitate the description of the magnetic structures 
(see below) we have redefined the triclinic basis vector 
c so that c — —2a' — c', where a' and c' are the basis 



vectors used in Refll3l. The important crystallographic 
features for the description of the magnetic properties 
are easily identified. The Fe'^+ ions, all in the high- 
spin state S=5/2 (see below), are arranged in clusters, 
separated by (V04)'^~ groups, containing non-magnetic 
ions. Each cluster of 6 Fe'^+ ions consists of two 
identical FesOia monomers, related by a center of inver- 
sion. With our convention, the c- axis runs along the 
line connecting adjacent clusters through their centers 
and open ends. In Fi^ we have also represented the 
two main types of magnetic interactions: intra-clustcr 
direct-exchange and SE interactions (thick lines), which 
are antiferromagnetic (AF), and inter-cluster SSE inter- 
actions (thin lines). There is clearly the potential for 
magnetic frustration, since Fig[T] reveals that the effec- 
tive inter-cluster and intra-cluster Fe-Fe linkages form 
loops (label 1 and 2), which contain an odd number of 
Fe sites, within which coUinear AF arrangements cannot 
be fully satificd. 

The magnetic frustration and competition between SE 
and SSE interactions give rise to a complex low tem- 
perature phase diagram, reminiscent of that of other 
cycloidal magnets. Above 100 K the magnetic suscep- 
tibility follows a Curie- Weiss law (not shown). The 
effective paramagnetic moment fJ-eff = 6.103(4)/iB is 
near the value expected for Fe'^+ ions in the high-spin 
S=5/2 state (At*^^ = 5.91^b). The Weiss temperature 
&CW = —124.9{i)K confirms the presence of strong AF 
interactions. The low temperature data clearly indicate 
magnetic transitions at Tjvi=22K and Tjv2=15K, consis- 
tent with a frustration index of \0cw\/Tni ~ 6. The spe- 
cific heat, (Fig. [5] a) also shows two lambda-type anoma- 



lies at these critical temperatures. The magnetic contri- 
bution to the specific heat was determined by subtracting 
the lattice contribution, indicated by the dashed line in 
Fig. [Ha), which was estimated by fitting the data at high 
temperature to the Debye function (Debye temperature 
of 385(5)K)). The magnetic entropy integrated between 
2K and 50K is Sm =13.98 J.K^^.mol^^, approaching the 
classical value of i?log(2S' + 1) =14.89 J.K-^.mol"^ for 
S=5/2. A large fraction of the total magnetic entropy 

30%) is only recovered far above T^vi, indicative of 
short-range AF ordering above T jvi [ill ■ 

To obtain further insight into the nature of the mag- 
netic ordering taking place at Tjvi and Tpf2 , magnetic 
measurements were made on the needle-shaped single 
crystals. The crystal orientation was set at 300K before 
cooling the crystal in a 3-T field and recording magneti- 
zation data between 5 and 35K. At low temperatures, the 
direction of minimal susceptibility coincides with the di- 
rection of growth of the needle-shaped crystals (approxi- 
mately along the crystallographic direction a). With this 
field orientation, labeled H|| in Fig. [^b, one only observes 
the magnetic transition at Tjvi. The transition at Tn2 
is clearly visible in measurements with the field applied 
perpendicular to the a axis, i.e., in the b* — c* plane. By 
rotating the crystal around its needle axis, we have de- 
termined the secondary easy magnetic directions in this 
plane. This is shown in Fig. The minimum and max- 
imum of the magnetization were found for the crystal ro- 
tated (j}min = —40° and (jfmax = 50° away from the c* di- 
rection, respectively. The temperature dependance of the 
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FIG. 2: Magneto-electric phase diagram for FeV04. (a) Heat 
capacity (filled circles). The dashed line shows the estimated 
lattice Debye contribution (see text for details) (b) Single 
crystal magnetization at 3 Tesla with the magnetic field ap- 
plied in the directions Hy (filled circles) and Hx (opened cir- 
cles) explicated in the panel (e). c) Bulk electric polarization 
measured for a sample cooled in a positive (E-|-,red) and neg- 
ative (E-,blue) electric field of 160V and dielectric constant 
(black solid line) d) Angular dependence of the magnetization 
at 3 Tesla and T=5K measured varying the orientation of the 
applied magnetic field in the (b*, c*) plane. H is aligned with 
c* when 6 = 0. 
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FIG. 3: (Color) Zoom in on the magnetic order of the Fe 5- 
shaped cluster of FigU containing the surrounded spin. The 
three inequivalent Fe sites are shown as small red, green and 
blue spheres, and the Fe sites with primed labels are ob- 
tained by inversion symmetry. The Fe are arranged in two 
types of edge-sharing polyhedra: two Fe'^'''05 trigonal bypi- 
ramids (sites 2-2') and the 4 other sites in Fe^^Oe octahe- 
dra. Spin are aligned with a common direction A at T=18K, 
and rotate in the same (A,B) plane at T=2K. The mag- 
netic moment amplitudes (black labels) were calculated using 
equation [T] with the refined parameters given in the supple- 
mentary information [itI] . accounting for a de-phasing term of 
68.4° coming from the position of the cluster in the crystal 
(R„ = (0 2 1)). Only the values of the refined phases tpj are 
given in degrees (red labels) for clarity. 

magnetization with the field orientation (pmin , labeled as 
in Fig. ^Bp shows a pronounced drop below TAf^flit. 
The phase boundaries indicated in Fig. [5] therefore mark 
the domains of stability of two distinct magnetic phases 
(I and II). The temperature dependence of the dielectric 
susceptibility, shown in Fig. [IJc), has no anomaly at T^ri 
but the sharp peak at Tjv2 and the concomitant appear- 
ance of an electrical polarization (Fig. 2c) demonstrates 
that the system becomes ferroelectric in phase II, while 
phase I is not ferroelectric. The value of the electrical po- 
larization at low temperatures is an order of magnitude 
smaller than for TbMnOsQ, a large reduction which can 
be partly attributed to the fact that our dielectric mea- 
surements were made on polycrystalline samples. 

The magnetic structures of phases I and II were 
determined by single-crystal neutron diffraction. Be- 
low Tati, the data show new Bragg peaks of mag- 
netic origin, which can all be indexed with a single, 
nearly temperature-independent propagation vector k = 
(0.222 - 0.089 0.012) almost perpendicular to c. In 
the most general case of a single-k incommensurately 
modulated structure, the magnetic moments Mj(Ri) on 
a given crystallographic site describe an ellipse as they 
propagate in different unit cells R;. Mj(R;) can there- 
fore be written as: 

Mj(Ri) = Aj cos(27r.k • Ri + (/3j) + Bj sin(27r.k • Ri + (/3j) 

(1) 

where Aj and B-,- are two perpendicular vectors defin- 
ing the major and minor semiaxes of the ellipse (isj. An 




FIG. 4; (Color) Magnetic structure of FeV04 as derived from 
refinements of the single crystal neutron diffraction data T — 
2K (top) and T = 18K (bottom). The magnetic moments 
are shown as brown arrows. For the non-coUinear structure at 
2K, the envelop of the helicoidal modulation is also shown by 
grey ellipses at each lattice site. Frustrated exchange loops are 
reproduced form Fig[T] and one spin on loop 1 experiencing 
the geometrical exchange frustration is surrounded. 



initial set of refinements indicated that the magnetic 
structure of phase I is a collinear spin-density wave (i.e., 
Bj= 0), whereas phase II possesses a helical structure. 
It is noteworthy that the collinear direction of the spins 
in phase I coincides with the major semiaxis Aj of the 
ellipses in phase II, indicating that the spin anisotropy 
does not change at the phase boundary. In principle, 
both the amplitudes and directions of Aj and Bj as well 
as the phase angles of the modulations can be different for 
each site. We have shown that the correlation between 
parameters can be significantly reduced by introducing 
constraints [l3] whilst still yielding refinements of excel- 
lent quality. This has lead us to consider only minimal 
models for phase I (10 parameters, Rf2 = 10.3% and 
= 1.7) and phase II (14 parameters, Rf2 = 5.07% 
and and = 8.2). The constraints retained suggest 
that the data do not support significantly different ori- 
entations of the ellipses and that the amplitudes | Aj | and 
Bjl (figlS]) of atoms related by centro-symmetry do not 
differ significantly, although their phases are not related 
by centro-symmetry and need to be determined indepen- 
dently. Perspective views of the magnetic structures of 
phases I and II are shown for one cluster in Fig. [3] and 
for few clusters in FiglH 

The helical magnetic order reveals different degree of 
frustration in different directions. It is characterized by 
the presence of quasi-lD AF order on chains of S-shaped 
clusters running in the c direction. The SSE path linking 
clusters in this direction is therefore probably the least 
frustrated SSE interaction inducing only very slow rota- 
tions of the average AF direction over a very long period 
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of approximately 110 nm. This contrasts with all the 
other SSE pathes, which induce large rotations of the av- 
erage AF direction between neighboring chains. 

Noteworthily, none of the magnetic orders is cen- 
trosymmetric because the difference between the phases 
of sites related by inversion symmetry deviates signifi- 
cantly from 180° and similar phases are refined at T=2K 
and T=18K (Figl3|). This inversion symmetry breaking 
is clearly visible at T=2K. Each half cluster of 3 Fe shows 
an almost perfect AF order, but between the two halves, 
the AF directions arc canted by 15 ~ 20°. The collinear 
magnetic order at T=18K simply appears as the projec- 
tion of the helical order on the common A direction. We 
argue that the observed k vector and the refined phases 
are here fixed by the symmetric part of the Heisenberg 
magnetic exchange energy independently of the presence 
or not of a B component, and that other energy terms 
control the stability of the collinear state with respect to 
the non collinear state. 

The transition to collinear magnetic order in multi- 
ferroic magnets with a spiral ground state has already 
been discussed by Mostovoy[10j. At the temperature 
T7V2, there is a cross-over between a high temperature 
regime in which anisotropy dominates so that the mo- 
ments align along the local easy axis, and one in which 
competition between anisotropy and entropy favors dis- 
ordered components of moment in directions orthogonal 
to it. For FeV04 the local easy axis is the A direction 
deduced from the neutron data, which coincides with the 
macroscopic easy axis, labelled Hy determined from the 
magnetisation measurements and labelled H II in Fig. 2b. 
In phase II the moments are confined to the plane of the 
ellipses with A and B as semi-axes, drawn in Figs 3 and 
4. When a field is applied in any direction orthogonal to 
A (H||) the temperature dependence of the magnetisa- 
tion shows a singularity at Tn2, a rapid fall below Tn2 
and almost no singularity at Tni- The orthogonal di- 
rection for which the fall is most pronounced is parallel 
to B, the minor axis of the ellipse. The non-coUinearity 
of phase II is due to the stabilization of the second order 
parameter modulating moments oriented parallel to B; it 
reflects the presence of additional terms in the free energy 
which outweigh the entropy effect. The non-coUinearity 
suggests that these terms are due to antisymmetric spin- 
orbit coupling of the Dzyaloshinski-Morya type propor- 
tional to vector products of spins S; x Sj . Most theories 
of multiferroic spiral magnets have invoked this interac- 
tion as the driving force leading to magnetically induced 
ferro-electricity. However it is not yet clear whether it 
is spin currents associated with these terms or the 
simultaneous emergence of ion displacements linked to 
configurational asymmetry in the magnetic structure, [l^ 
which give the larger contribution to the ferroelectric po- 
larisation. 

In summary, we have demonstrated magnetically in- 
duced ferroelectricity in FeV04. It is noteworthy that its 



magnetoelectric phase diagram in the vicinity of Tn2 is 
very similar to that of the classical TbMnOa system in 
the vicinity of T/ocfe=27K[3, [i^l- This gives a strong in- 
dication that spin-orbit interactions may make an impor- 
tant contribution to the electrical polarization in FeV04. 
However, the Fe^~^ ion has L=0 in the free-ion state and 
in FeV04 should be orbitally quenched which reinforces 
the induced nature of the mechanisms responsible for 
multiferroic behavior in spiral magnets. The magneto- 
electric coupling is simpler in FeV04, because the only 
magnetic atoms participating are Fe; it therefore not in- 
fluenced or complicated by rare earth magnetism, as this 
is the case in TbMnOajloi. 

This study of FeV04 has clarified the relative roles of 
exchange frustration, anisotropy, entropy, and antisym- 
metric exchange in defining the magnetic state of "type- 
11" multiferroics. In TbMnOs it is difhcult to identify 
the origin of the magnetic frustration because the SSE 
interactions form a network of next nearest neighbor in- 
teractions interpenetrating that of the first neighbor su- 
perexchange interactions. In FeV04 on the other we have 
established a clear connection between the moment re- 
duction or the spin rotations, which are typical of the 
incommensurate magnetic phases of spiral magnets, and 
the most relevant SSE paths giving rise to frustration. 

We would like to thank Dr J. A. Rodriguez- Velazman 
for his support during the neutron diffraction experiment. 
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In the case of FeV04, the modulation and the symmetry of magnetic structures is more conveniently described in 
an alternative cell setting (a, b,c) — (a',b', —2a' — c'), where (a',b',c') was the cell used in the original structural 
description given in Refjlj. In the original work, cell parameter at RT are a' = 6.71633 A, b' = 8.06817 A, c' — 
'o6'' 9.34715 A, a' = 96.6537°,/3' = 106.528°, 7' = 101.499°, so that with our new setting, a = 6.71633 1, b = 8.06817 A, 
O ; c = 14.01312 A, a = 74.432°, (3 = 140.248°, 7 = 101.499°. The transformed Fe atomic positions in this new ceU are 
• given in table HI 

O 

Q 



Oh! 



> 



X 



Atom X y 



Fei -0.06433 0.69426 -0.40839 

Fe2 0.04169 0.89000 -0.21254 

Fes 0.05455 0.69432 0.01143 

Fe4 -0.05450 0.30573 -0.01143 

Fes -0.04169 0.11000 0.21250 

Ph ■ Fee 0.06443 0.30584 0.40839 
O 

o 

I ■ TABLE I: Atomic positions for FeV04 in the transformed cell (a, b, c) (see text) 



^ I The magnetic moments Mj(Ri) of any site j in a unit-cell translated by a vector Ri with respect to the origin is 
written: 

Mj(Ri) = Aj cos(2.7rk • Ri + ipj) + Bj sin(2.7rk • Ri + tp^) (1) 



' O ' where the constraint (Aj±Bj) is always imposed in our refinements. In FeV04, the propagation vector is k = 
^ (0.222 - 0.089 0.012). 

2 In the program FuUprofQ, we can constrain the direction of the parameters describing the magnetic structure. 

This is best done using spherical coordinates attached to a {x,y, z) orthogonal frame and defined in a way such as x 
coincides with a and z coincides with c*. 

For amplitude modulated phases (Bj=0) (table HH) . we use the conventional latitude and azimuth spherical angles 
'j^Oj) with ((j)j,9j) — (0,0) giving z, and {(f)j,9j) ~ (0,90°) giving x, in order to parameterize and/or constrain the 



Aj directions . 

. For helicoidal structures, we used the option in FuUprof which orients the vectors (Aj,Bj) and the vector Nj normal 

' to the (Aj,Bj) plane along the axis of orthogonal coordinate systems (x', y' , z') obtained by applying active rotations 

\ given by the Euler angles {(j)j,6j,Xj)i which rotates the (x,y,z) coordinate frame in turn along z, y, and z again. 

I The parameters Xj are in fact always constrained to be the same for all sites to avoid correlations with the phases ipj 

00 ' (table Ell). 

O 



^ ■ [1] B. Robertson and E. Kostiner, J. Solid State Chem. 4, 29 (1972). 
5^ i [2] J. Rodriguez-Carvajal, Physica B 192, 55 (1993). 



Model I: fixed directions A 
for all sites and symmetry constraints on phases 
(Rb = 38.9% = 32.4 : not converging) 
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Model III: Model I + 




relaxation of the synmietry constraints on phases 
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FCr, 


1-4.1 


<t>r, = (hi 


Or, = Oi 


-17(6) 






W(j = '.'i 


(h, = Oi 


-201(2) 


Model IV: model 1 + 



relaxation of the spin directions (partial) 
and phases symmetry constraints 





[Rb 


= 10.2% 


= 1-62) 




Atom 


|A| 


0(°) 


0n 


o 


Fei 


3.00(5) 


3(2) 


95(2) 





Fe2 


3.23(5) 


7(2) 


94(2) 


-168(9) 


Fe3 


2.88(4) 


1(2) 


97(2) 


3(2) 


Fe4 


1^3 1 <, 


i>4 = 03 


04=03 


-197(2) 


Fes 


\M <, 


= 02 


05 = 02 


-20(9) 


Fee 


\M <, 


i>6 = 01 


06 = 01 


-203(2) 



TABLE II: Magnetic parameters refined at T=18K 



3 



Model I: fixed directions for all sites and symmetry constraints on phases 
{Rb = 27.9% = 270.0: not converging) 



Atom 


1^1 \B\ 


xC) 


0(°) 








FEi 


4.0(5) 4.6(5) 


XI = -55(5) 


= 87(6) di 


= 41.3(6) 





FE2 


3.1(5) 3.3(5) 


X2 = Xl <? 


1)2 = 01 


^2 = Ol 


188(3) 


FEs 


5.0(4) 3.3(5) 


X3 = Xi </ 


^3 = (61 


Oa = Ol 


21(2) 


FE4 


\M \B3\ 


X4 = Xl <f 


!)4 = 01 


O4 = Ol 


ip:i + 180 


FE5 


\M \B2\ 


X5 = XI <f 


is = 01 


05=01 


ifi2 + 180 


FEe 




X6 = Xl <? 


!>6 = 01 


Oe = Ol 


ifii + 180 


Model II: model I + relaxation of the phases symmetry constraints 






{Rb = 5.07% 


X' = 8.24) 






Atom 


1^1 \B\ 


x(°) 


0(°) 


0{°) 




FEi 


4.51(7) 3.67(8) 


Xl = -77(2) 


= 81(2) 01 


= 33.2(1) 





FE2 


4.18(6) 3.16(7) 


X2 = Xl 9 


!)2 = 01 


O2 =0i 


-171(2) 


FE3 


4.29(7) 3.39(9) 


X3 = Xl <? 


!>3 = 01 


O3 = Ol 


7(7) 


FE4 


|A,| \B3\ 


X4 = Xl <f 


i4 = 01 


04 = Ol 


-189(7) 


FEs 




X5 = Xl <? 


is = 01 


05= Ol 


-18(2) 


FEe 


1^1 1 \Bi\ 


X6 = Xl <? 


!>6 = 01 


Os = Ol 


-195.2(7) 


Model III: Model II + partial relaxation of the directions constraints 




{Rb = 


6.29% x^ = 12.7 


: slightly diverging) 




Atom 


\A\ \B\ 


x(°) 


0(°) 


0{°) 




FEi 


4.52(8) 3.7(1) 


XI = -76(2) 


76(2) 


20.7(5) 





FE2 


4.21(8) 3.06(9) 


X2 = Xl 


73(2) 


28.8(3) 


-167(9) 


FE3 


4.31(8) 3.4(2) 


X3 = Xl 


85(3) 


42.6(3) 


0(7) 


FE4 


\A3\ \B:i\ 


X4 = Xl <? 


U = 03 


O4 = O3 


-196(7) 


FEs 


\A2\ IB2 


X5 = Xl <? 


is = 02 


Or, = O2 


-15(9) 


FEe 


l^il \Bi\ 


X6 = Xl <? 


ie = 01 


Oe = Ol 


-196.2(9) 



TABLE III: Magnetic parameters refined at T=2K 



